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SUMMARY 


A one-dimensional compression system model diat simulates the 
unsteady aerodynamics in the NASA J85-13/P^ (Planar Pressure Pulse 
Generator) test installation has been developed. In addition, the frequency 
response characteristics of the test apparatus were evaluated using a 
linearized frequency analysis technique. This analysis was performed in 
support of pelecting insti imentation locations to avoid standing wave nodes 
of the installation and thus, low signal levels. The feasibility of employing 
a set of explicit analytical surge criteria in evaluating the stability of an 
axial compressor was also investigated. 

The J85-13/P^ simulation program was developed by modifying an 
existing program known as the Dynamic Digital Blade Row Compressor 
System Stability Model (also called A^ROSTAP). The capability of the model 
to handle unsteady flows was enhanced by introducing an analytical expression 
for the unsteady-blade -lift force to coripute the blade force acting upon the 
fluid. The validity of the modified model was verified by comparing the 
simulation results with the test data for the FlOl Fan/P^G - a test installa- 
tion similar to the J85-13/P^G. 

Subsequent to the program validation, the model was run at 100% 
compressor design speed over the P^G frequency range of 100 to 850 Hz to 
simulate the stationary operations of the J85-13/P^G. Pressure waveforms 
at various locations of the test installation were obtained to provide a direct 
comparison with future test data. Then, the model was throttled to surge at 
five different frequencies, i« e. 100, 150, 200, 250, and 300 Hz. and the 
inlet planar wave distortion sensiHvity of the J85-13 turbojet engine was 
determined as a function of the P '^G frequency. 

The frequency response analysis of the J85-13/P^G was performed for 
the frequencies ranging from 10 to 1000 Hz using a General Electric develop- 
ed linearized frequency analysis program. This analysis yielded the nodes 
and anti- nodes of the test installation as a function of frequency. The 
resonance frequencies of the test apparatus were also identified from die 
analysis. 

The explicit analytical surge criteria were derived from the linearized 
conservation equations (mass and momentum) based upon the Lyapunov 
first theorem of stability. The criteria were applied to two different con- 
figurations, ''Moss'' and "Mehalic", of the J85-13 engine. They were found 
eff.*ctive at 100% desiga speed, but were ineffective at lower speeds. This 
is believed mainly due to the boundary conditions whose validity is weakened 
as the speed is lowered. Recommendations to improve the capability of the 
criteria are listed in the last section of this report. 
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2.0 INTRODUCTION 


NASA Lewis Research Center is developing a test installation known 
as die J85-13 /p3g, Figure 1, for the purpose of investigating unsteady aero- 
dynamics in gas turbine engines. The apparatus includes a General Electric 
J85-13 turbojet engine and a device known as the Planar Pressure Pulse 
Generator (P^). The P^* which was modelled after a similar device used 
earlier by General Electric (Reference 1), was designed and fabricated by 
Battelle Columbus Laboratories. 

The function of the P^G is to generate sinusoidally oscillating flow over 
a wide range of frequencies and amplitudes. This is accomplished by 
modulating the flow area with a rotating disc and a stationary disc with mating 
holes in each. Hie frequency of the planar wave is varied by changing die 
rotor speed. The amplitude of die oscillation is varied by changii^ the 
rotor -to -stator spacing. When die J85-13 engine is operated coupled with 
the P^G, the engine can be subjected to oscillating inlet ilows in which both 
die frequency and amplitude are independently variable over wide ranges. 

The J85-13/P^G test is the first of a series of experimental efforts 
planned at NASA Lewis Research Center. Hie objective of this test is to 
evaluate the unsteady aerodynamic responses of a multi-stage compressor* 
e. g. unsteady-operating -point excursion, compressor instability, etc., 
with oscillating inlet flow conditions. The J85-13 ei^ine was selected for 
this test because of the extensive operational and analytical experience diat 
NASA has obtained with the compressor. A large aunount of data has been 
accumulated for both the performance and distortion sensitivity of the J85-13 
compressor (References 2 and 3). In addition, several computer models 
have been developed to simulate die unsteady as well as the steady operations 
of the compressor (References 2, 3, and 4). 

The main objective of the effort discussed in this report was to develop 
a digital computer simulation code that is capable of modelling the unsteady 
behavior of a compressor subjected to the oscillating inlet flow generated 
by the P^G. This program will be utilized for the pretest predictions as 
well as for posttest analyses to aid in obtaining an understanding of die 
unsteady responses of the compressor. In addition, the frequency response 
characteristics of the J85-13/P^C test facility were studied in support of 
selecting instrument locations to avoid the standing wave nodes within die 
apparatus and thus low signal levels. The feasibility of employing explicit 
analytical expressions for surge prediction was also investigated as part of 
this effort. A set of explicit surge criteria was developed as a result of this 
study. These criteria are applicable only «.o an isolated rotor although the 
inlet and exit rotor boundary conditions are obtained from the simulation of a 
complete compression system including stators and free volumes. 
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The J85ol3/P^ simulation pr<^ram was developed by modifying the 
Dynamic Digital Blade Row Compression System Stability Model (called 
AEROSTAP, Reference 5) developed by die General Electric Company. This 
program employs a one-dimensional, pitch line, blade row-by-blade row 
model to solve die diree conservation equations widi the proper thermo- 
dynamic relationships. These diree governing equations (mass, momentum, 
and energy) are written in a volume -average form and the dynamic solution 
is obtained through a time-marching scheme. Derivation of the three govem- 
it^ equations and the solution method employed in the AEROSTAP Program 
is outlined in Sub' ction 3. 1. 'Die essence of die modifications to this pro- 
gram is the incorporation of an expression for the unsteady-blade lift force 
that is induced by the oscillating inlet flow produced by the P^. The 
expression employed to describe the unsteady-lift force of a cascade was 
derived from the unsteady-lift force diat i& applicable to an isolated thin 
airfoil subjected to gust loading (Reference 6). The derivation of this 
expression is presented in Subsection 3. 3. 

The validity of employing the resultant expression for the unsteady- 
blade -lift force for a cascade was verified by comparing the simulation 
results with experimental data obtained during a test similar to that proposed 
for the J85-13/P^G (Reference 7). This verification procedure and the 
res’ilts are discussed in Subsection 4. 1. The modified AEROSTAP Program 
was then utilized to simulate the operations of die JSS-IS/P^G at various 
operating conditions. The results of the simulation are presented also in 
Subsection 4. 1. 

The frequency analysis of the J85-13/P^G test facility was performed 
using a digital frequenc, analysis program known as the DSP (Digital Surge 
Prediction) program (Reference 4). In this program, the frequency response 
characteristics are evaluated by linearizing the three conservation equations 
and then determining the transfer function of the linearized system by means 
of Laplace transforms. These techniques which are discussed in detail in 
Reference 4 are outlined in Subsection 3. 4. 

The analysis of the frequency response characteristics of the J85-13/ 
P^G was conducted by dividing the geometry into two regions; i.e. , the 
region forward of the P^G and the region aft of the P~G including die J85-13 
engine. This analysis yielded the resonance frequencies of the two sections 
and their higher-order harmonics as well as die nodal points at various 
frequencies. These results are presented in Subsection 4. 2. 

The explicit analytical stability criteria were derived based upon the 
Lyapunov's first theorem; the dieorem states that a system of nonlinear 
equations is stable if the eigenvalues of the linearized version of the system 
have negative real parts. In this derivation, the eigenvalues of the linearized 
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conservation equations are examined by tiie Routh>Hurwits criterion. At 
this point of development and for die ease of analysis, only the mass and 
momentum conservation equations are employed; i. e. the flow dirough a 
compression system is assumed homentropic. The explicit stability 
criteria are derived in Subsection 3. 5. A more detailed derivation of the 
equations is given in Appendix A. Application of diese criteria and a 
critique of the technique are found in Subsection 4. 3. 
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3. 0 ANALYTICAL TECHNIQUES 


ITiis section presents the analytical techniques used to develop the 
digital computer program that simulates the operation of the J85«13/P^ 
test installation. A brief description of the techniques employed to invest!* 
gate the frequency response characteristics of die J85-13/P^ apparatus is 
also presented. In addition, derivation of explicit analytical stability criteria 
for axial compressors is discussed. 

The J85*13/P^ simulation program was developed by modifying a 
version of the Dynamic Digital Blade Row Compression System Stability Model 
(called AEROSTAP, Reference 5) developed by the General Electric Company 
to simulate a test facility, known as the FlOl Pan/P^, similar to that of the 
J85*13/P^G. Reference 2 indicated diat while this program adequately 
simulates most aspects of FlOl Fan/P^G operations, the phase angles of the 
pressure waveforms computed at various locations in the test facility did not 
compare well with the test data. It was believed that this mismatch in phase 
angles is due to inadequate description of the unsteady behavior of the com* 
pressor blades under the oscillating inlet flow generated by the Planar 
Pressure Pulse Generator. Thus, a novel method of describing the unsteady 
behavior of the compressor cascade has been added to the AEROSTAP to 
enhance the accuracy of the simulation. 

Previously, the axial force acting on the fluid due to the compressor 
blades was computed via the tangential blade force which was obtained from 
the generalized Euler Turbine Equation. This approach has been utilized 
successfully for the simulation of steady as well as transient operations of 
various compressors. However, the expression for the tangential force does 
not permit accounting for the phase difference between the oscillating flow 
incident on the blade and the similarily oscillating blade force. Therefore, in 
the new approach, the axial blade force is computed through the blade-lift 
force within which the phase difference can be easily included. 

The governing equations of motion that are in the form used in the 
AEROSTAP are derived in Subsection 3. 1. The ancillary relations required 
to solve the governing equations and the solution technique are also presented 
in the subsection. Modelling of the P^G which generates sinusoidally 
oscillating flow over a wide range of frequencies and amplitude is discussed 
in Subsection 3. 2. 

The expression for the unsteady-blade-Uft force is derived in Sub- 
section 3. 3. This expression was derived from the unsteady-lift force that 
is applicable to an isolated thin airfoil subjected to gust loading. It was 
assumed that by introducing correction constants, the resultant relationship 
is valid also for a cascade. 
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Ilia fraquancy analysis of ^a J85->13/P^G tast facility was conduetad 
with the ttsa of a frequency analysis program known as tiia Digital Surge 
Prediction (DSP) program. The frequency response characteristics are 
evaluated by linearizing the tiiree conservation eqtiations and subsequently 
determining the transfer function of tiie linearized system by means of 
Laplace transforms. This program is also capable of determining the 
stability of an equilibritun solution of the three conservation equations by 
examining the eigenvalues of tiie linearized system. NASA developed sub- 
routines are used in &e DSP program for the evaluation of Laplace '.rans- 
form and for die inspection of die eigenvalues. This program which is 
discussed in detail in Reference 4 is briefly described in Subsection 3.4. 

In Subsection 3. 5. the explicit analytical criteria are derived based 
upon die Poincare » Lyapunov stability dieory. These criteria check die 
stability of an equilibrium solution of the diree governing equations by 
examining the eigenvalues of the linearized conservation equations by the 
Roudi-Hurwitz criteria. These stability criteria are applied to each rotor 
volume of a compressor. It is assumed that the operation of the compressor 
is unstable if the sign of the real part of die eigenvalues associated with one 
or more of the rotor volumes changed from negative to positive. 


3. 1 ANALYTICAL TECHNIQUE FOR THE SIMULATION OF COMPRESSION 
SYSTEM 


In this subsection, the equations of change which describe the motion 
of the compressible fluid flowing through a compression system are derived. 
They are written in volume-averaged form in terms of die nomenclature 
given in Pages xiii through xiv. In addition, the auxiliary relationships used 
in the consarvation equations are also derived. 

The three conservation equations, mass, momentum, and energy, are 
solved by a technique known as the time-marching scheme. This technique 
which is discussed in detail in Reference 3 is outlined in Subsection 3. 1. 2. 

3. 1. 1 Governing Equations 

The equations of motion are derived for an arbitrary volume, V|^, with 
the inlet area, A|, the exit area, Aj 4 . i, and the volume lengdi, shown 

in Figure 2. The macroscopic balances of mass, momentum, and energy 
for this volume for quasi one-dimensional flow yield the following equations: 
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Continttity; 




*itil 


( 1 ) 


Momentum: 




[ 

w? 1 



L^Vi 

*c‘’l*l*i+I * * «*> **' * Kj 


Energy: 


dPS 

at 



Vi 





(J> 



Figure 2. Quasi One>Dimensional Volume. 

where the quantities with tiie bar " - " are defined by the following relation- 
ship. 




( 4 ) 
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This set of die conservation equations provides an exact description of 
the quasi one -dimensional flow dirough an arbitrary volume. Note that the 
only assumption applied in the derivation of these equations is diat the flow is 
one-dimensionaL All the boundary effects such as the wall friction loss and 
heat transfer can be included in the right-hand side of die momentum and 
energy equation dirough £Fic and Sp. However, it should be noted that diese 
boundary effects are not included in the solution since their orders of 
magnitude are much smaller than odier terms in the conservation equations. 

In order to obtain a solution of die diree conservation equations, it is 
necessary to supply the ideal gas equation of state, the caloric equations of 
state and expre. ns for SFk end Sp. 

The force ...ms, in Elquation 2 represents the force exerted on 

the fluid by all the solid surfaces such as the end walls and blades. As stated 
previously, the wall friction is neglected in diis analysis. Then, die force 
term can be expressed as 

XF,^= /pdA + F^ (5) 

lateral 
area 


where dA is the differential area vector and F is the axial-direction force 
acting on the fluid due to the compressor blades. 


The lateral surface pressure integral can be evaluated by 

fpdA = fpdA 

lateral 1 
area 


( 6 ) 


Introduction of the definition of a mean pressure 



enables one to rewrite £>)uation 6 as 


/pdA ^ 
*lateral 
area 


"m<^i - 


(7) 


( 8 ) 
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Three different expressions for Pj^ are given in References 1 and 4; 1) a 
polynomial representation developed from tiie isentropic relations for blade* 
free volumes 


•■m ■ I “• (Mj - M,) + »• 000908125 

2 X 

- 0.24098565 - Mj) + 0.02419741 (mJ-mJ) 

+ 0. 03442452 (M® - M ] (9) 


2) an empirical relationship for blade*free volvunes 




where C = Cj(A 2 - Aj)/Aj + 


Cj, = empirical constants 


and 3) an empirical correlation for bladed volumes 


■ 




( 10 ) 

( 11 ) 


( 12 ) 


where Pjj is the higher of the values of the inlet or exit static pressure and 
Pl is the lower of the two. The empirical correlation for blade -free 
volumes was developed by examining the behavior of the mean pressure in 
the blade-free volumes of the J85-13 engine. The correlation for bladed 
volximes was derived from the observation of P^ in bladed volumes of 
various General Electric compressors where there were no losses. This 
empirical relationship is assumed to be equally valid when blade losses are 
present. 


The other term in the sum of the surface forces is the axial force, F,, 

£ 

exerted by the blade. This force can be deternuned from the various force 
components depicted in the following sketch. 
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Figure 3. Components of Blade Forces. 

Examination of die force vectors in the sketch shows diat 
F^sF^tan^. . Fj^^ 


where Ff is the tangential force at the pitch line and Fq^ axial 

component of the blade drag force acting upon the fluid. The lift direction 
angle* in the expression of F 2 is given by the following equation. 


p, » <e, 

(14) 

where + 1 

(15) 

(Sz = » 

(16) 


: lift correction angle 

/3*, fi*-. metal angles at the inlet and exit respectively 


The tangential force, F<p, can be obtained from the generalized Euler 
Turbine Equation and is written as 



J ^ 


(17) 
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where r j and C 2 are the pitch line radii at the inlet and exit of the blade row, 
respectively. The second term in Equation 17 represents the tangential 
force due to the storage of angular momentum within a blade row volume. 

The drag force acts perpendicular to the lift force at an angle of 
from the axial direction. Thus, the axial component of the drag force, 
can- be written as 


*■0. * 


where 



1 


(18) 


(19) 


The prime (') symbol indicates that the quantities are evaluated w*th respect 
to the relative velocity frame of reference. The area, , is the flow area 
perpendicular to the direction of the inlet relative air angle, j3j. The 
relative total -pressure loss coefficient, < 0 ^ is defined by the equation. 


at - 


^2) 


t2 'Ideal 


•Kz) 


Actual 


P • P 
tl 1 


( 20 ) 


The values of are evaluated from the experimental data that are correlated 
as a polynomial function of tan i. 


The term Sp in Equation 3 is known as the internal entropy generation 
term in the absence of the heat transfer at the boundaries* This term can be 
expressed as 


WR In 
F 


) Actual 


( 21 ) 


The ideal total-pressure ratio can be written (Reference 8 ) as 



^tl /ideal 




( 22 ) 


12 



(23) 


where * 2 wNr 2 /a^j 

I 

and < 4 ^ is the inlet speed of sound evaluated with the relative total 
temperature. The expression for the actual total-pressure ratio can be 
obtained by rearrangement of Equation 20 and substitution of die isentropic 
relation for total pressure and static pressure as 



For more detailed derivatior of the governing equations and t.ieir 
ancillary relationships, readers sh Id consult Reference 3. 

3.2.2 Solution Technique 

The technique employed in solving the governing equations. Equations 1 
through 3, is known as the time marching scheme. This technique is based 
upon a Taylor series expansion of the volume -averaged flow properties at 
time, t, to predict the flow properties at time, t + At. Hie time derivatives 
of volume-averaged properties that are used in the Taylor series are 
obtained from the properties at the inlet or exit station of a vol.une. A 
suitable interpolation scheme should be selected to assure the numerical 
stability of the solution technique. Reference 3 gives a detailed account and 
a critique of various interpolation schemes; thus, they are not discussed here. 

Once the volume-averaged properties and their time derivatives are 
o' cained, a Taylor expansion of the variables, P, W, and pS is straight 
f.x'ward. This expansion is illustrated for one variable, namely, the 
volume -averaged density. The Taylor series for P correct to second order 
can be written as 

P(t + At) s P(t) + At + (25) 

dr ^ 

where P(t) is established by the initial conditions or a vious time step 
and 
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( 27 ) 



Examination of Equation 27 reveals that the right-hand side is composed of 
derivatives of station properties, and These derivatives are 
evaluated from the momentum equation through an interpolation scheme since 
Equation 2 provides only die volume -aver aged derivative. Substitution of 
Equations 26 and 27 into Equation 25 permits one to evaluate the volume- 
averaged density correct to second order at the next time increment. _ 
Similarily, diic technique can be used for the remaining two variables, W and 
7S, and can be continued from one time step to the next for the desired 
number of time steps. 

This computational technique is illustrated in Figure 4 in block diagram 
format. For more detailed discussion of the solution technique, readers 
should consult Reference 3. 


3, 2 Modelling of Planar Pressure Pulse Generator 

The function of the Planar Pressure Pulse Generator (P^G) is to 
generate a sinusoidally oscillating flow over a wide range of frequencies and 
amplitudes. This is accomplished by modulating the flow area with a 
rotating disc and a stationary disc with mating holes in each. The air flow is 
induced by the compressor downstream of the P^G. The frequency ot the 
planar wave is varied by changing the rotor speed while the amplitude of the 
oscillation is varied by changing the rotor-to- stator spacing. The change in 
the rotor-to-stator spacing alters the effective flow area and, thus, varies 
the amplitude of the oscillation. 

The P^G area is essentially choked at all times and the total pressure 
loss results mainly from shock losses as the supersonic stream immediately 
after the P^G adjusts to higher pressure further downstream. Some viscous 
losses are present due to mixing of the air streams from adjacent holes in the 
stator. However, for the flew rates and planar wave amplitude 'f interest, 
the shock losses are predominent; therefore, the viscous losses ^ > neglected 
in this model. 


The P^G model approximates the aerodynamics of the actual operation 
as follows: A steady-state condition satisfying the mean conditions of the 
P^G operation is established using the mean steady-state total-pressure 
recovery across the P^G as meastired during the test. This enables one to 
compute an effective choked area at the P^G from the flow function 
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FF » 



where A is the choked area. The Clow Conction and ti^e total temperature 
across die remains constant. This mean area becomes the base for die 
amplitude of the area fluctuation. The amplitude is controlled by the para» 
meter "amp" as illustrated in die following equation. 

^3G ' ^3g base [ » ♦ (2»Nt/60) ] (29) 

This choked area can be related to die upstream duct A* /A and provides 

the exit boundary condition to the P^G. 

As much of die existing General Electric Dynamic Digital Blade Row 
Compression Component Stability Model logic was utilized as possible, and 
the P^ model was constructed as a subroutine to the main program. For 
further information concerning the modelling and operation of die the 
readers should consult References 1 and 7. 

3. 3 Blade Lift Force 

In Subsection 3. 1. the axial-direction blade force, F 2 , was expressed 
in terms of the tangential forces, F -p, and the axial component of the drag 
force, FDz* stated previously, this approach of computing Fj was highly 

successful in simulating transient as well as steady-state operations of axial 
compressors under various operating conditions. However, when a large 
amplitude oscillation such as the planar pressure wave generated by the 
P^G ic present at the inlet, this method cannot satisfactorily describe the 
frequency response, namely the phase angle, of the compressor (Reference 7). 
Therefore, a new method of computing F^ has been added to die AEROSTAP 
program to improve the frequency response of the simulation. This new 
technique calculates the axial-directicn blade force through the blade lift 
force, Fp,, instead of the tangential force, Fx; the drag force, Fd, is 
evaluated in the same manner as described before. From the force compon- 
ent diagram depicted in Figure 3, the new relationship can be derived as 

Fz s Fl siu/Sa - F£)COs3» (30) 

The expression for the unsteady-blade-lift force, Fp,* was derived 
from die lift force applicable to a thin airfoil subjected to sinusoidal vertical 
gust. That is: First, the unsteady lift force for the thin airfoil is obtained 
by considering the lift force due only to the gust loading. It is assumed diat 
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the airfoil is placed initially in a uniform incompressible velocity field of U 
without any angle of attack. Next, the uniform velocity field is tilted slightly 
to give a steady-state lift force to file airfoil. This is intended to approxi- 
mate file flow over a compressor blade. The total lift force acting on fiie 
fluid by the airfoil is fiien obtained by superposing the unsteady lift force 
determined previously on fiie steady lift force. Subsequently, this expression 
for the total lift force is applied to a cascade by introducing two correction 
factors. These factors account for fiie cascade effect, profile shape and 
camber, profile drag and compressibility. 

For a fiiin airfoil. Figure 5, encountering a sinusoidally oscillating 
vertical gust of Wq = WQei^^ in an otherwise uniform stream of U, 

Reference 6 expresses the unsteady lift force due only to the gust loading as 

^'l. dyn =WPUbW^ |c (K) | Jo(K) - | + (K) } e^"“* (31) 



Figure 5. Thin Airfoil Under Sinusoidal Vertical 
Gust, Initial Angle of Attack is Zero. 

where fiie Theodorsen function, C (K), is defined by 

C(K) - F(K) + jG(K) 




( 32 ) 
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Hq (n : 0, 1) in right-hand side cf Equation 32 is the Hankel function of 
tile second kind which is a combination of Bessel functions of tiie first and 
second kinds: 



(K) = J (K) + jY^(K) 
n n 


( 33 ) 


The complex argument^ K, is given as 


K = 


ftib . lib 
2U ■ ^ 2U 


(34) 


For constant amplitude oscillations, however, |tis zero and K is expressed 
simply as 


K = 



(35) 


Sub jtittttion of tiie Theodorsen function into Equation 31 and rearrange- 
ment gives 


F, , = irPUbW_ (A + jB)e^***' (36) 

Lig uyn O 

where 

A = F(K)Jo(K) + G(K)Ji(K) (37) 

B = Ji(K) +G(K)Jo(K) - F(K)Ji(K) (38) 


This expression for the unsteady lift force can be rewritten with the aid of a 
trigometric relation as 


DYN ej*^*e^^, tf=tan'^B/A) (39) 

Examination of Figure 5 reveals that if Wq«Up the time -dependent 
angle of attack, aQg can be expressed as 

jait 

“G = ® (40) 
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Sobstitation of tiiis angular relationship into E^iuation 39 yields 


*■ 1 , dyn * + («) 

Note tiiat ^e bracketed term in ^e right-hand side of the above expression 
is equivalent to the steady-state lift force except diat the angle of attack, Oq, 
varies sinusoidally widi time. The term, -t- accounts for dte change 
in the amplitude of &e lift force due to tiie transient effect. The {diase 
difference between the gust velocity fluctuation and tiie lift force oscillation 
is represented by the \mit vector, e ifi . 

A similar derivation can be applied to the case where a steady-state 
lift force exists prior to the introduction of gust loading. The steady-state 
lift force is induced by giving the airfoil an initial angle of attack, org^. 
Figure 6. 



Figure 6. Thin Airfoil Under Sinusoidal Vertical 
Gust, Initial Angle of Attack is 

Under the flow conditions illustrated in the above sketch, the total lift force 
can be obtained by superposing the unsteady lift force due to a gust on the 
steady-state lift force. Reference 6, i. e. 


F 


L 


L, ss 


+ F 


L, dyn 


where 


F, = vpU^ba 
L, ss ss 


(42) 


(43) 


It is assTuned that Ogg is small. Since ocq is also assumed to be small, the 
dynamic lift force represented by Equation 41 is equally valid for this case. 
Thus, the total lift force, F^, becomes 
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A* + 


Now, ^is erpression is applied to a compressor blade by introducing 
two correction factors, K j and K 2 , and replacing Oq by tiie incidence angle 
fluctuation, I = Ai e That is 


F, s bK,« 

L 2 ss 


[‘*K^ 


+ B^e 


-1 *^ 1 ® 

where 9 = tan ) 

A 

As stated previously, Kj ax^ provide corrections for the cascade effect, 
profile shape and camber, profile drag and compressibility. 

The steady^state lift forces 

F, s sAU^bK^o (46) 

L» S8 2 S8 

for the compressor blades are computed in the AEROSTAP program witii the 
proper K 2 for each blade. The steady- state^lift force for each blade at a 
given corrected speed was computed and represented as a polynomial 
function of tan i where i is the incidence angle. In order to utilize diese 
polynomials* Ek^uation 45 is rearranged as 


-r f 

L L, 8s[ 


I A^ 4 b" Ale^ 


L, ss 


+ d)j 


This expression can be rewritten by substituting the relationship 
e^^ = cos <6 + j sin<^ where ^ ~ o»t+ $ 


''l. ss I ‘ - Va^ + B^ Ai[cos («t +d) + j sin(eit + « )] 
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Ekiuation 49 represents the final expression for the unsteady»blade-lift 
force. Note that the dynamic portion of &e lift force is composed of real and 
imaginary parts. Use of die real or imaginary part in die actual computation 
depends upon die oscillatory mode of 1 = STeJ^^. If die incidence angle 
fluctuates as a sine wave, the imaginary term is used widiout j; while real 
part is utilized if I changes as a cosine wave. 

Evaluation of A and B in Equation 49 widi a digital computer is difficult 
since Hankel and Bessel functions are not readily available as a library 
functions. Consequently, die functions, F (K) and G (K), Figures 7 and 8, as 
well as the first and second order Bessel functions of the first kind are 
curve-fitted by polynomials and included in AEROSTAP as a subroutine. The 
amplitude of die incidence angle fluctuation, STI is computed from diat of 
die mass flow oscillation produced by the P^. 

Finally, it should be noted that the correction factor, K| was evaluated 
by comparing the simulation results with die test data of the FlOl Fan/P^ 
installation. As will be shown later, Kj was found to be very close to unity. 

3.4 FREQUENCY RESPONSE ANALYSIS TECHNIQUE 


The frequency analysis technique employed in the Digital Surge Pre- 
diction (DSP) program is essentially a Laplace transfer analysis of a set of 
linear equations. 


dx. 

1 


dt 


= A..X. + B. 
ij 1 in 


U 

n 


R = C .X. 
m mi 1 

where i = 1, 2, 3, 
j = 1. 2, 3. 
n = 1, 2, 3, 
m= 1, 2, 3, 


. 3K 
. 3K 
. L 

. P 


K: Number of model volumes 
L: Number of forcing function variables 
P: Number of output variables 


(50) 

(51) 
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Equation 50 describes the sensitivity of tiie perturbations, X{, of the 
independent variables about an equilibrium solution to ^e forcing function. 
Bin ^n* 1^8 equation is obtained by linearising the three conservation 
equations written in vector form as 


dVi 

-3T = ' 


(52) 


and adding the forcing function, Bi,^ Uq. Note that the conservation equations 
are applied to all the volumes md die volume-averaged density, mass flow 
rate and entropy for all the volumes are taken as independent variables. 
Therefore, diere are 3K equations for 3K independent variables. Lineari- 
zation of Equation 52 about an equilibrium solution, y^j, yields 

dy. 

cvi-v.,) (5J) 

} 

Introduction of a matrix 


A.. 

ij 



(54) 


and the definition of the perturbation about y 


X. 

1 


y. - y . 
1 ei 


into Equation 53 gives 


dx. 

= A .. X. 
dt ij 1 


(55) 


(56) 


Now, addition of the forcing function, Uq, to this relationship yields 
Equation 50. The forcing function is formulated from the inlet and exit 
boundary conditions or any other conditions imposed within the model. 

Equation 51 represents the output response in terms of a linear com- 
bination of the perturbations. The vector, is the output variables such 
as the pressures, mass flow rates, etc. at various locations in the model. 
The matrix, expresses the linear relationship between the output 

variables and ^e perturbation variables. 
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The Laplace transform equivalents of die system represented by 


Equations 50 and 51 ara 



(57) 


(58) 


where s is the Laplace variable. These two equations can be combined and 
rearranged in transfer hinction form as 


1 

it ' 7 V = C .(si.. - A..) *B. 

U (s) mi ij ij in 


(59) 


where I., is the unity matrix, 
ij 

Equation 59 forms the basis for the frequency response analysis. Eval- 
uation of this relationship with a digital computer is. however, not desirable 
since it involves the inversion of a matrix that contains a free parameter. 
Therefore, each element of the transfer function is evaluated separately in 
the DSP program. That is, the transfer function between the K-th element 

of R and 1 -th element of U 
m n 

l{Rj^| det (W..) 

rTu7r det (sl.^ - A..) 

* 1* ij ij 

for the desired combination of m and n are computed individually. 

mat *ix, W is a Wronskian Matrix whose derivation is explained in detail in 

Reference 4. 

The determinants in the transfer function represented by Equation 60 
are each polynomials in s. Therefore, the frequency response is numerically 
evaluated by substitution s = jci> into these polynomials where 6»is a given 
frequency and j is the imaginary number, 

3. 5 EXPUCIT AN^.LYTICAL SURGE CRITERIA 


One of the important parameters in die design and operation of gas tur- 
bine engines is surge margin. If a gas turbine engine is operated in the 
compressor stall region, not only is the engine performance adversely 
affected, but also the engine life deteriorates due to high mechanical 
stresses and high temperatures. Thus, the normal operating line is selected 
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below the surge line with proper allowances for tiie external as well as the 
internal destabilizing effects such as inlet distortion, power transients, 
thernal acceleration transient, manufacturing tolerances, etc. However, 
lowering of the operating point from ^e surge limit usually results in 
decreasii^ the pressure ratio and efficiency of the engine. Therefore, the 
necessity of determining the exact stall line and subsequently, the proper 
surge margin fiiat assures efficient and successful steady or transient 
engine operation is evident. 

Stall lines for compressors have been determined usually through 
ei^;ine tests. However, this process is costly and time consuming. Conse- 
quently, there is a need to obtain stall lines analytically. In recognition of 
this need, fiie General Electric Company has developed two independent 
techniques for estimating compressor stability limits. These two analytical 
methods are discussed briefly below. 

The first method utilizes the time -dependent solution of the General 
Electric developed Aerodynamic Stability Analysis Program (AEROSTAP, 
see Subsection 3.2). This method determines the stability limits by 
examining two time -dependent parameters (Reference 3), i. e. , 1) The 
amplification factor, (dW/dt) rotor ^^*^^^***^exit mass flow ratio, 

(W2/W1) rotor- 


The first parameter represents the ratio between the time-rate change 
of mass flow rate of a rotor and that of the exit volume of the model. It was 
observed that the amplification factor becomes greater than two in all the 
rotor volumes at the experimentally determined stability limit. The second 
parameter is the ratio of the exit mass flow rate to the inlet mass flow rate 
of a rotor volume. At the experimental surge limits, this ratio for one or 
more rotors deviates from unity by more than 0. 00004. This means that a 
compressor system becomes unstable if more than 0. 004% of the inlet flow is 
stored or depleted from any of the rotor volumes. 

These two stability criteria based on the time -dependent solution of the 
AEROSTAP simultaneously indicate instability at the experimental surge 
limits. Therefore, they are used in support of each other. 

The second method of predicting stall line utilizes the steady-state 
solution or a quasi -steady- state equilibrium solution of the AEROSTAP. This 
technique is based on the Poincare - Lyapunov theory of stability. This 
theory states that an equilibrium solution of any nonlinear system of the form. 

V ^n> = 2,....., n) (61) 
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Based on this principle, die three conservation equations express in 
volume-averaged form (see Subsection 3.2) are linearired to obtain the 

square matrix, (i> K = 1, 2 n), ^.here n is equal to 3 x 

number of volumes. Then the eigo'*' values of this matrix are examined by 
die Routh-Hurwitz criterion to determine the stability of the solution of the 
conservation equations. One or more of the eigenvalues will have a positive 
real part at the stability limit. 

The two methods described thus far have been utilized successhilly in 
computing the stability limits of various General Electric manufactured gas 
turbir. ' ei^ines (References 3, 4, and 7). The results show an excellent 
agreement between the computed stability limits and thr experimentally 
determined surge lines. 

At the present time, these two mediods provide adequate estimates of 
the surge limits. However, these techniques are tedious and costly; the 
first method requires a slow quasi-steady-state throttling along a speed line 
and the latter technique reqrires evaluation of a large size matrix and its 
characteristic equation. Therefore, there is a need to develop a new tech- 
nique that is simple and less expensive. Further, it is desirable to possess 
a technique that can readily show the relationships between the occurrence of 
the instability and the aerodynamic and geor etric characteristics of a com- 
pressor. In recognition of diese needs, explicit analytical stability criteria 
for the evaluation of the surge limits were derived. 
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Tliis new techni<{ue, detcribed hereinafter, is also derived from tiie 
Poincare « Lyapunov tiieory of stability. However, this metiiod diHers signi* 
ficantly from the previously described method in the fact tiist die st^ility is 
accessed for individual rotors; whereas die previous mediod evaluates the 
stability of a compressor as a whole. FOrdier, for the ease of analysis, the 
flow tiirough the compression system is assumed homentropic. 


For die homentropic flow of a compressible fluid in a rotor volume. 
Figure 9, the governing equations can be written, in volume>averi^ed form, 
as 


dt 


( Py W^) 



(64) 


dW 

dt 



l»l. W^) 



1 



‘cVl 




+ PjAj - PjAj 


- P, 


M 


(Ai . A^) + Fg 




where P, . = 
M 


— P + P 
3 1 3 2 


( 65 ) 

( 66 ) 


The boundary conditions are 

P^j, = const, (at inlet) 


( 67 ) 


W 



= const, (at exit) 


( 68 ) 


These boundary conditions are based on the assumption that any small per- 
turbations in the rotor volume under consideration cannot propagate into the 
adjacent rotor volumes due to an abrupt change in acoustic impedence 
between rotors. This assumption is quite valid at high speeds: however, its 
validity is weakened as the operating speed is lowered. 
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Derivation of the partial derivative8« Equations 70 ttiroi^h 73. dF^/dP, 
and is given in Appendix A. * 


According to the Poincare •Lyapunov stability dieory. die solution for 
the simultaneous differential equations. Equations 64 and 65. are stable if 
all die eigenvalues of die matrix. Equation 69. have negative real parts. The 
characteristic equation of diis matrix is 
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-»s 0 


(74) 


The Roudi'Hurwitz criterion states diat die roots. X» of this equation have 
negative real parts if all die coefficients of this polynomial are positive. 
Therefore, the stability criteria are expressed as 


ap. 


af. 


aw. 


< 0 


af. af. 
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aw. 


af. 


aP, 


(75) 

(76) 


These two inequalities are the final expressions for the stability criteria 
and can be evaluated widi die substitution of partial derivatives given in 
Appendix A. Application of these criteria to a stability evaluation of the 
J6S-13 "Mehalic" and "Moss" engines and die results are discussed in Sub- 
section 4. 3. A critique of die method is also presented. 


31 



4.0 PRESENTATION OF RESULTS 


This section presents tiie results of tlie analysis performed using tiie 
analytical techniques outlined in foe previous section. First, the results of 
the J8S*13/P^ simulations are discussed in Subsection 4. 1. This subsection 
includes a description of foe digital model of foe J85>13/P^ test apparatus, 
validation of foe simulation program, time-dependent analysis of foe test 
installation and an estimate of foe stability limit of foe J85 compressor wifo 
inlet planar wave distortion. Second, foe frequency analysis of foe test 
apparatus is described in Subsection 4. 2, Items discussed in fois subsection 
are foe amplitude transmission and foe nodal analysis of foe system model. 
Finally, foe results obtained by applying foe explicit analytical stability 
criteria to two different configurations of foe J85 compressor are presented 
in Subsection 4. 3 along with a critique of the criteria. 

4.1 J85-13/p3g simulation 

The model geometry, foe procedure for operating the model, and 
verification of foe model capabilities are discussed in the following para- 
graphs as are the results of time -dependent simulation of the J85-13/P^ 
operation and surge line estimations. 

4. 1. 1 Geometry and Operating Procedures of foe J85-13/P^G Model 

The geometry of foe J85-13/P^G test installation simulated by the 
digital model was shown previously in Figure 1. The entire setup from the 
bellmouth entrance to the combustor exit is modelled by 79 volumes. Table 1 
defines the individual characteristics of these 79 volumes. The area distri- 
bution of the test apparatus as used in this simulation is shown in Figures 10 
and 11. Note that the volume representing the plenum (Volume 1) and the 
first five volumes of the bellmouth (Volumes 2 through 6) were replaced by 
smaller constant area volumes, denoted by the dotted line in Figure 10. 

This was necessary since the large areas as well as the large area gradients 
in this region caused numerical instability in the AEROSTAP program. 

The physical characteristics of the compressor geometry, i. e. , metal 
angles, pitch-line radii, areas and solidities, are defined in Table 1 for each 
blade row. 

The locations of the instrumentation are shown in Table 1. Instru- 
mentation locations that will be used during the J85-13/P^G test are 
duplicated in the model to provide a direct comparison between the model 
results and the future test data. The parameters chosen for this comparison 
are defined in Table 2. 
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Table 1. System Geometry 




Station 



Area (In.^) 


Vobune Description 

Number 

Instrumentation 

Physical 

Blockage 

EffecH 

Plennm 

1 


16286.00 

1.0 

16286. 

Bellmoath Lip 

2 


2827.43 

1.0 

2827. 

1 

1 

3 


2135. so 

1.0 

2135.. 



4 


1851.42 

1.0 

1851. 



5 


1659.74 

1.0 

1659. 



6 


1397.08 

1.0 

1397. 

Bellmoutb 








7 


1219. 22 

1.0 

1219. 



8 


1090.73 

l.O 

1090. 



9 


995.72 

1.0 

995. 



10 


925. 52 

1.0 

925. 



11 


874. 90 

l.O 

874. 



12 

PTBME, PSBME 

840. 64 

1.0 

840. 



13 


820. 82 

1.0 

820. 

Belli 

nouth EIxit. Plane 0« 0 

14 


814. 33 

1.0 

814. 



IS 


814. 33 

1.0 

814. 



16 

PSFPI 

814. 33 

1.0 

814 



17 


814.33 

1.0 

814 


framb 


tv, ' •* - -i. 


|py)OUT 


le 1. Svrtem Geometry 


Physical 

Area (In.^) 
Blockage 

Effective 

RR (b.) 

Vohune 

Number 



Length 

Solidity 

16286. 00 

1.0 

16286. 00 

-- 

1 

— 


4.3 


2827.43 

1.0 

2827. 43 


2 



1.7 


2135. 50 

1.0 

2135. 50 

— 

3 

— 


2.0 

m m 

1851.42 

1.0 

1851. 42 

— 

4 

-- 

— 

2.0 

• • 

1659. 74 

1.0 

1659. 74 

— 

5 


-- 

4.0 

•- 

1397.08 

1.0 

1397.08 


6 

-- 


4.0 


1219. 22 

1.0 

1219. 22 

— 

7 


-- 

4.0 

— 

1090.73 

1.0 

1090.73 

— 

8 

-- 

-• 

4.0 

-- 

995.72 

1.0 

995. 72 

— 

9 

— 

-- 

4.0 

• - 

925. 52 

1.0 

925. 52 

— 

10 

-• 

-- 

4.0 

-- 

874. 90 

1.0 

874. 90 

— 

11 

M * 

-- 

4.0 


840, 64 

1.0 

840. 64 

-- 

12 

-- 

-- 

4.0 

-- 

820. 82 

1.0 

820. 82 

— 

13 

-- 


4.0 

-• 

814.33 

1.0 

814. 33 


14 

-- 

-- 

4. 24 


814. 33 

1.0 

814. 33 

— 

15 

— 


4.0 

• • 

814.33 

1.0 

814.33 

— 

16 

A A 

» mi 

4.0 

-- 

814. 33 

1.0 

814.33 

— 

17 

-- 

mm m 

2. 56 
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Volume Description 


Duct Ahead of P^G 


Duct Ahead of P^G 




Station 

Number 

Instramentation 

Physical 

Area (In.^) 
Blockage 

Effectiv 

18 


839. 82 

1.0 

839. 8: 

19 

PSFP2 

897. 27 

1.0 

897. 2‘ 

20 


-951.15 

1.0 

951. 1. 

21 


1006. 60 

1.0 

1006.6 

22 


1063. 62 

1.0 

1063.6 

23 

PSFP3 

1086. 86 

1.0 

1086.8 

24 


1052.65 

1.0 

1052. 6 

25 


1001.92 

1.0 

1001. : 

26 


991.83 

1.0 

1 1 

• 

27 

PSFP4 

991.83 

1.0 

991.: 

28 


991.83 

1.0 

991. 

29 

PSFP5 

966. 27 

1.0 

966.2 

30 


856.11 

1.0 

856. 

31 


689. 23 

1.0 

689... 

32 


506.68 

1.0 

506. 

33 


404. 51 

1.0 

404., 

34 


315.09 

1.0 

315. 

35 


263. 15 

1.0 

263. 


iKAMB ! 

PRECEDING PAGE BLANK NOT FILMED 


Physical 

Area (In.^) 
Blockage 

Effective 

RR (In. ) 

Volume 

Number 

(Deg. ) 

(Deg.) 

Length 

fei 

Solidity 

839. 82 

1.0 

839. 82 

— 

18 

-- 

^ * 

4.0 


897. 27 

1.0 

897. 27 


19 

-- 

— 

4.0 


951.15 

1.0 

951.15 

» A 

20 

— 

— 

4.0 


1006. 60 

1.0 

1006. 60 

— ^ 

21 

— 

-- 

4.0 


1063. 62 

1.0 

1063.62 

— 

22 

— 


2. 55 


1086. 86 

1.0 

1086. 86 


23 

— 

— 

4.0 


1052.65 

1.0 

1052.65 

— 

24 

A ^ 


4.0 


1001.92 

1.0 

1001.92 

-- 

25 

« * 

— 

4.0 

* ^ 

991.83 

1.0 

991.83 


26 

•- 

-• 

4.0 


991.83 

1.0 

991.83 

— 

27 

— 

— 

4.0 


991.83 

1.0 

991.83 

— 

28 

— 

— 

4.4 


966. 27 

1.0 

966. 27 

— 

29 

— 

— 

2.0 


856. 11 

1.0 

856. 11 

— 

30 

— 

A ^ 

2.0 


689. 23 

1.0 

689. 23 

— 

31 

A » 

— 

2.0 


506.68 

1.0 

506. 68 

— 

32 

— 

— 

2.0 


404. 51 

1.0 

404. 51 

-• 

33 


— 

2.0 


315.09 

1.0 

315.09 

-- 

34 

-- 

-• 

2.0 

• • 

263. 15 

1.0 

263. 15 

— 

35 


-- 

2.0 



35 

> 

C 







Station 



Area (In. 


Volume Description 

ft 

Number 

bistrumentation 

Physical 

Blockage 

Effec 



36 


229. 30 

1.0 

22 

Duct Ahead of P^G 








37 


210.20 

1.0 

21 



38 


209.51 

1.0 

20« 



39 


204.92 

1.0 

20^ 



40 


203.83 

1.0 

2<: 

P^G Rotor Inlet Plane 








41 


198. 67 

1.0 

i: 




PSFP6, PTFP6 




Discharge Plane 

42 


192.54 

1.0 

1 










43 

PSAPI, PTAPI 

192. 54 

1.0 

1' 



44 


192. 54 

1.0 

1 



45 


192. 54 

1.0 

1 



46 


192. 54 

1.0 

1 

Duct Ahead of J85 Compressor 







47 


192. 54 

1.0 

1 



48 


206.46 

1.0 




49 


187. 96 

1.0 


PlatM 

s 2.0 

SO 

PSAP2 

187.96 

1.0 



^ i 








51 


187. 96 

1.0 




52 


187. 96 

1.0 


Engii 

le Inlet IGV 

53 

PTIGV, WHGV 

187. 11 

1.0 



PRECEDING PAGE BLANK NOT RLMED 



Physical 

Area (In. 
Blockage 

Effective 

RR (In.) 

Volume 

Nnmber 

(DegJ 

(p«»-i 

Lengtii 

(In.) 

Solidity 

229. 30 

1.0 

229. 30 

-- 

36 



2.0 

A • 

210. 20 

1.0 

210.20 


37 

-- 

« m 

2.0 

W « 

209. 51 

1.0 

209. 51 

— 

38 

— 

•• 

2.0 

mm w 

204.92 

1.0 

204. 92 

-• 

39 

— 

-- 

2.0 


203. 83 

1.0 

203. 83 

— 

40 

-- 


2.0 

-• 

198. 67 

1.0 

198. 67 

— 

41 



2.75 

-- 

192. 54 

1.0 

192. 54 

— 

42 

-- 

-- 

1.0 

— 

192. 54 

1.0 

192. 54 

— 

43 

-- 


2.0 


192. 54 

1.0 

192. 54 

-- 

44 


-- 

2.0 


192. 54 

1.0 

192. 54 


45 

-- 

-- 

2.0 

-• 

192. 54 

1.0 

192. 54 

— 

46 

-• 

-- 

2.0 

-• 

192. 54 

1.0 

192. 54 

— 

47 

• • 

-- 

2.0 

-- 

206.46 

1.0 

206. 46 

•• 

48 

•- 

-- 

2.0 


187.96 

1.0 

187. 96 

-- 

49 

-• 


1.0 

•• 

187. 96 

1.0 

187. 96 

— 

50 

-• 

* 

1.0 

-- 

187. 96 

1.0 

187. 96 

-- 

51 

-- 

*• 

2.0 

-• 

187. 96 

1.0 

187. 96 


52 

-• 

M m 

2.25 


187. 11 

1.0 

187.11 

5.918 

53 

OeO 


2.51 

3. 000 
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fOUX)UT FRAMB 



Volume DftBcriPtiott 


Free Volume 
R1 

51 
R2 

52 
R3 

53 
R4 

54 
R5 

55 
R6 

56 
R7 

57 
R8 

58 

OGV 


Station 

Number Instrumentation 

Physical 

Area (In.^) 
Blockaee 

Effe< 

54 

165. 53 

0.97 

16 

r 

55 

153.48 

0.97 

1 

56 

136. 38 

0.97 

13 

57 

120. 89 

0.97 

11 

58 

108.48 

0.97 

1 

59 

98.23 

0.97 


60 

90.42 

0.97 


61 

81.07 

0.97 


62 

75.07 

0.97 


63 

67. 54 

0.97 


64 

63. 45 

0.97 


65 

58.40 

0.97 


66 

56. 04 

0.97 


67 

53.25 

0.97 


68 

52.47 

0.97 


69 

51.55 

0.97 


70 

51. 01 

0.97 


71 

51. 01 

0. 97 



PReC£D.NG PAGE BUNK NOT 

rotooux I 


Physical 

Area (In.^) 
Blockage 

Effective 

Volume 
RR (In.) Number 

(S*sJ 

(Deg. ) 

Length 

(»»•) 

Solidity 

165. 53 

0.97 

160. 56 

6. 195 

54 


0.0 

0.92 

- • 

153.48 

0.97 

148. 88 

6.272 

55 

48.39 

37.5 

1.186 

1.143 

136. 38 

0.97 

132. 29 

6.400 

56 

34. 84 

5.4 

1.041 

1.057 

120. 89 

0.97 

117. 26 

6. 504 

57 

50. 20 

30. 20 

0. 789 

1.204 

108.48 

0.97 

105. 23 

6.600 

58 

40. 58 

12. 10 

0.848 

1.208 

98.23 

0.97 

95. 28 

6,712 

59 

51.05 

26.49 

0. 599 

1.326 

90.42 

0.97 

87.71 

6.792 

60 

45.68 

14.80 

0.740 

1.350 

81.07 

0.97 

78.64 

6.896 

61 

50.81 

26.69 

0.517 

1.384 

75. 07 

0.97 

72.82 

6. 960 

62 

46. 40 

19.40 

0.655 

1.390 

67.54 

0.97 

65.51 

7.032 

63 

55.95 

28. 35 

0.430 

1.432 

63.45 

0.97 

61.55 

7.080 

64 

46.03 

23. 20 

0.617 

1.436 

58. 40 

0.97 

56.65 

7.128 

65 

54.80 

34.10 

0. 363 

1.306 

56. 04 

0. 97 

54. 36 

7. 152 

66 

44. 90 

22.20 

0. 584 

1.311 

53.25 

0. 97 

51.65 

7. 184 

67 

55.11 

40. 37 

0. 333 

1.318 

52.47 

0.97 

50.90 

7.192 

68 

43.79 

21.10 

0. 597 

1.322 

51.55 

0.97 

50.01 

7. 200 

69 

52.94 

44.82 

0. 303 

1.326 

51.01 

0.97 

49.48 

7. 208 

70 

42. 10 

9. 58 

0. 466 

1.209 

51.01 

0. 97 

49.48 

7. 205 

71 

9. 58 

0. 00 

0. 417 

1.209 
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Station 



Area (In.^) 


Volume Description 

Number 

Instrumentation 

Physical 

Blockage 

Effect! 

Compressor Discharge 

72 

PSACl 

51.01 

0. 97 

49. 



73 


65. 19 

1.0 

65. 

Combustor Inlet 

k 

74 


124.66 

1.0 

124. 

1 

1 

75 


182. 53 

1.0 

182.. 



76 


189. 10 

1.0 

189.. 

Combustor Volume 








77 


189. 10 

1.0 

189. 



78 


184.83 

1.0 

184. 



79 


94.87 

1.0 

94. 

Turbine Nozzle Inlet 

80 


97.95 

1,0 

97. 


^CEDING PAGE 


blank Nor n 




iiUOiA 


Physical 

Area (In.^) 
Blockage 

Effective 

RR (In.) 

Volume 

Number 


(P^sj 

Length 

(to.) 

Solidity 

51.01 

0.97 

49.48 

7.205 

72 



1.0 


65. 19 

1.0 

65. 19 

• ^ 

73 

— 

— 

2.0 

• «» 

124. 66 

l.O 

124. 66 

— 

74 


— 

2.0 

• • 

182. 53 

1.0 

182. 53 


75 

— 

— 

2.0 


o 

• 

00 

1.0 

189. 10 


76 

-- 

— 

2.0 


189. 10 

1.0 

189. 10 


77 

-- 


2.0 


184.83 

1.0 

184. 83 

•- 

78 

— 

— 

1.5 


94. 87 

1.0 

94. 87 


79 

-- 


1.0 


97.95 

1.0 

97. 95 
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Table 2. J85-13/P^G Instrumentatton Location 


Nanvt 

Location 

Reference Plane 

PTBME 

1 

• 

q 

Plane 0. 0 

PS3ME 

1 

• 

q 

Plane 0.0 

PSFPl 

+ 12.1'' 

Plane 0. 0 

PSFP2 

+ 21.6" 

Plane 0.0 

PSFP3 

+ 35. 5" 

Plane 0. 0 

PSFP4 

+ 51.7" 

Plane 0. 0 

PSFP5 

+ 59. 5" 

Plane 0. 0 

PSFP6 

+ 82.8" 

Plane 0. 0 

PTFP6 

+ 82.8" 

Plane 0.0 

PSAPl 

- 14.0" 

Plane 2. 0 

PTAPl 

1 

• 

q 

Plane 2. 0 

PSAP2 

0.0" 

Plane 2. 0 

PTIGV 

5.3” 

Plane 2. 0 

wncv 

+ 5.3" 

Plane 2. 0 

PSAC! 

+ 19.2" 

Plane 2. 0 

PR-OA 

Overfall Pressure Ratio 

a) First two characters PT, PS and WI stand for total pressure. 

static pressure and physical mass Pow rate. 

respectively. 

b) « and 
plane, 

+ denote upstream and downstream from the refcxence 
, respectively. 

Table 3. FI 01 Fan/P^G Instrumentation Location 


Name 

Description 

P^GPS 

Static pressure at the P^G strut 

AXP53 

Static pressure at the fan inlet 

PTIGV 

Total pressure at the inlet guide vane 

PS14 

Static pressure at the fan exit 

FDPSC 

Static pressure at the splitter 
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The performance characteristics of the eight stage compressor are 
represented in the form of blade row relative total-pressure loss coefficients. 

and deviation angles, as a function of incidence angle, tan i. In these 
representations, all losses are assigned to the rotor and it is assumed that 
only lossless turning of flow occurs in the stators. The loss coefficients and 
deviation angles of the J85 compressor are discussed in detail in Reference 3. 

The blade characteristics. «ui'and 8, represent blade-row performance 
up to the steady-state surge line. However, due to the large amplitude 
planar waves that are to be imposed on the compression system, the incidence 
angles encountered during model operation may instantaneously exceed the 
steady-state surge line. Therefore, an extrapolation technique is employed 
to extend the blade characteristics into the post steady-state surge regime. 

A detailed discussion of this extrapolation is also given in Reference 3, 

The system model is subject to two boundary conditions; i. e. , constant 
total-pressure and total-temperature condition at the bellmouth inlet and 
constant fiow function at the combustor exit. In an actual engine test, the 
entire test setup or the bellmouth is placed in a large plenum. Usually, the 
total pressure and total temperature of the plenum are not affected by the 
operation of the engine. Therefore, in the digital model, a plenum volume 
is added in front of the bellmouth and the constant total-pressure and 
total -temperature condition is employed as the boundary condition. During 
the actual operation of a gas turbine engine, the turbine diaphragm is choked 
due to high pressure and high temperature of the combustor gas. This 
choked condition is simulated by assigning the constant flow function boundary 
condition to the combustor exit plane. 

The method of simulating the dynamic events of the J85-13/P^G with 
the model is very similar to the actual test procedure used for FlOl-Fan/P^ 
test. Reference 7. The test procedure was as follows: First, the compres- 
sion component w'as brought to steady- state operation at a low operating point 
on a chosen speed line. Then, the P^G was actuated with the rotor-to-stator 
(RTS) spacing at its maximum. This was done to maintain the minimum 
amplitude of the planar wave while the P^G accelerates to a desired rotational 
speed (frequency). After the P^ had attained the desired frequency, the RTS 
spacing was decreased to increase the planar wave amplitude to a desired 
value. Similarily, the model enters the computation at the chosen operating 
point on a speed line and simul||tes the steady operation at this point for a 
period of time. This process is/called "settling”. Then, a frequency is 
assigned to the P^G and the amplitude of the planar wave is increased 
gradually by increasing the amplitude of the choked area fluctuation (see 
Subsection 3. 2). This process is called "ramping". After the ramping, the 
model can be run with the chosen frequency and amplitude of the planar wave 
in a stationary manner and subsequently, throttled to an instability point. 
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After the J85>13/P^G model geometry axnl the model operation scheme 
were established, the newly modified AER OSTAP was run wifii the FI 01 Fan/ 
P^G geometry as well as wifii the J8S<*13/P^ geometry for the purpose of 
validating the computer program. These validation runs and file results are 
discussed in the following subsection. 

4.1.2 Program Verification 

Veiification of the program capabilities was accomplished in a threefold 
manner: 1) Evaluation of one of the correction factors used in the expression 
for fite unsteady-blade 'lift force, 2) Verification of the validity of employing 
the lift force in describing die aerodynamic behavior of compressor blades, 
and 3) Assessment of the numerical stability of the AEROSTAP program 
when it is run widi the J85-13/P^ geometry. 

As shown in Subsection 3. 3, the expression for the unsteady-blade-lift 
force contains two correction factors, Kj and K 2 . These factors correct the 
expression for the cascade effect, compressibility, profile shape and camber, 
profile drag, etc. The first correction factor, Kj, provides an adjustment to 
the phase angle of the dynamic portion of the lift force while K 2 provides 
corrections to the steady- state lift force as well as to the amplitude of the 
dynamic lift force. 

As stated in Subsection 3.3, K 2 is accounted for when the steady- state 
lift force is computed within the AEROSTAP program. This is accomplished 
through the lift direction correction angle, /3g, Reference 3. The values of 
fif. are determined experimentally from compressor performance data. 
Therefore, only Kj is evaluated in this study. 

The phase difference between the oscillating incidence angle and blade 
lift force was defined in Subsection 3. 3. 

$ = tan'^ (Ki A/B) (77) 

A and B in the above relationship are given in terms of the zeroth and first 
order Bessel functions of the first and second kinds and are plotted as a 
function of K = iub/ZXJ in Figures 12 and 13, respectively. The values of d for 
Kx = 1 are plotted in Figure 14 also as a function of K. 

In order to determine the correction factor, experimentally, it is 
necessary to obtain the phase angle, $, from test data. However, evaluation 
of 0 is extremely difficult, if not impossible, at this point since there is no 
reliable means of deducing the lift force from the FlOl Fan/P^G test data; 
however, the incidence angle can be computed readily from the experimental 
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data. Tlxerefore, evaluation of K| was attempted through a comparison of 
the pressure waveforms obtained during ^e FlOl Fan/P^ test with those 
simulated with K| as an independent variable. The experimental as well as 
the theoretical waveforms were obtained at five different locations in die 
FlOl Fan/P^ test installation (Figure 15). Tliese locations are defined in 
Table 3. 

Operation of the FlOl Fan/P^ wifii a P^ frequency of 118 Hz and an 
IGV total>pressure oscillation of 15% peak*to>peak about the mean pressure 
were simulated for the Kj values ranging from 0. 5 to 1. 5 in increments of 
0. 1. The average corrected mass flow rate and total-pressure ration across 
the fan were 151.8 Kg/sec (334.6 Ibjj^/sec) and 2.243. respectively. The 
pressure waveforms computed at the five different locations were compared 
with the corresponding experimental waveforms. The result indicated that 
the best agreement between the two sets of data was achieved when K| was 
between 0.9 and 1.1. These K| values yielded little differences in die 
simulation results. 

A similar comparison was performed at a P^G frequency of 350 Hz to 
determine the effect of frequency on K^. Other operating conditions were 
similar to those for the 118 Hz case. The result of this comparison also 
showed that the best agreement between the test data and simulation results 
was obtained when Kj is close to unity. 

In addition, the relationship of K| to the amplitude of the IGV total- 
pressure variation was investigated. The amplitudes of the IGV total- 
pressure oscillation for both the 118 Hz and 350 Hz runs were increased to 
22. 8% and 25.6% peak-to-peak about the mean pressure, respectively. The 
result of this investigation indicated that Kj is also independent of the 
amplitude of the inlet *otal-pressure fluctuation. 

Therefore, it was concluded that the correction factor, is a constant 
and is independent of the frequency and the amplitude of the planar inlet dis- 
tortion. The constant value of Kj was chosen to be unity since a small 
variation of K| about unity resulted in little difference in the simulation 
results. 

After the correction factor, Kj was chosen, operation of the FlOl Fan/ 
P^G was simulated at five different frequencies, i. e. , 42, 80, 118, 220, and 
350 Hz, for the purpose of verifying the validity of employing the lift force in 
describing the aerodynamic behavior of compressor blades. Previously, the 
AER OSTAP program computed the blade force acting on the fluid from the 
blade tangential force and the drag force. The newly modified AEROSTAP 
utilizes the blade lift force and the drag force in computing the blade force. 
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The validity of employing die new approach cannot be demonstrated easily 
through any theoretical deductions. Therefore, die validation was attempted 
through a comparative study of the experimental data widi simulation results 
at die five frequencies previously listed. 

In Figures 16 through 20, the results of the simulation are superposed 
on the experimental pressure waveforms at various locations in die test 
apparatus. The solid line represents the experimental data while the dotted 
line represents die simulation results. Examination of these figures 
indicates that the two sets of data are generally in a good agreement with 
each other within computational and experimental errors. Therefore, it was 
concluded that the use of the unsteady-blade-lift force in evaluating *^r.e blade 
force is a valid means of describing the aerodynamic behavior of a compres- 
sor blade. 

The remaining discrepancy between die experimental and theoretical 
waveforms may be attributable to the following computational factors. 

1) Model does not include the higher order harmonics of the 
P^G that are present in the test data (amplitude and phase 
angle). 

2) One -dimensional volume averaged model with a finite volume 
length (amplitude and phase angle). 

3) Model does not provide any dissipatory mechanisms such as 
viscous dissipation and heat transfer for the pressure waves 
(amplitude). 

Further, the polarity of the last two channels for 220 Hz run. Figure 19, 
might have been inverted during the recording or play back of the test data. 
The recording and play back were conducted with two different setups on 
different run days. 

Some of the results of the FlOl Fan/P^G simulations are presented 
graphically in Appendix B. 

Subsequent to the verification rxins with the FlOl Fan/P^G geometry, 
file AEROSTAP program was run with the J85-13/P^G geomet. y to assess 
the numerical stability of the program with the new geometry. 

At the beginning of this study, the model encountered a numerical in- 
stability that originated from reverse flow in the bellmouth region. Initially, 
the axial velocity in the bellmouth inlet region is very snaall, and the static 
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Figure 20. FlOl Fan/P^G Simulationf 100% 350Hz 



pressure is very close to the inlet total pressure, due to the extremely large 
flow area of this region. When die planar wave generated by the reaches 
this region, the flow reverses since the small kinetic head of the flow cannot 
overcome die pressure hill of die wave. Once die inlet velocity of a volume 
reverses, the density and die pressure of the volume decreases rapidly to 
satisfy the continuity relationship. This drop in pressure forces the velocity 
to the inlet of the adjacent downstream volume to reverse and subsequently, 
a similar process continues to the downstream volumes. Eventually, the 
reverse flow propagates down to the compressor and as a result, the compres- 
sor is "starved" and is thrown into an unstable operation that leads to the 
solution failure of the model. 

This type of instability, however, may not occur in actual test unless 
die amplitude of die planar wave propagated upstream to the bellmouth is 
unusually h*gh. If die area is too large for a given mass flow rate, the flow 
does not necessarily follow die contour of the bellmouth. Instead, an annular 
stagnation region exists near die solid wall and the air flows only through the 
core region of the bellmouth. Thus, the kinetic head of the flow can be high 
enough to counter the pressure rise due to the oscillation. 

In order to approximate the actual flow at the bellmouth, the plenum 
volume and the first five volumes of the bellmouth were replaced by smaller 
constar.t area volumes as described in the previous subsection. After this 
modiflcation to the geometry, the digital model was free of the instability due 
to reverse flow. 

With the modified geometry of the J85-13/P^G test aoparatus. the 
AEROSTAP program was run according to the model operation scheme 
established in the previous subsection. The model entered computation at the 
steady-state operating point at = 19.73 Kg /sec (43. 5 Ibjji/sec), see 
Figure 21. The total-pressure ratic at this point is 7.61. After "settling" 
for 1000 time steps, the P^G frequency of 100 Hz was assigned to the model 
and then the amplitude of the P^G area oscillation was "ramped" until the 
IGV total-pressure variation reached approximately 10% peak-to-psak about 
the mean total pressure. This ramping procedure is denoted by the dotted 
line in Figure 21. When the ramping was completed, the model was allowed 
to run steadily with a constant amplitude of the P^G area fluctuation for 
10, 000 time steps which are equivalent to 10 cycles of the ITO Hz oscillation. 

It was observed that the operating points of the J85 compressor formed a 
repetitive orbit of an ellipse. Figure 21, after three cycles. This clearly 
demonstrated that the model can be brought to a stable operation according 
to the operating procedure given in Subsection 4. 1. 1. 
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Subsequent to the program verification, the AEROSTAP prc^ram was 
run at eight different frequencies to investigate the time -dependent 
operation of fiie J85-13/P^ test installation and dien, fiirottled at four 
different frequencies to determine the instability points of the J85 compressor 
under the inlet planar distortions of various frequencies* The results of 
this series of investigations are presented in the following two subsections* 

4* 1. 3 Time-Dependent Analysis 

The time-dependent analysis of die J85-13/P^ test installation was 
performed at 100% compressor design speed at eight different frequencies, 
i.e., 100, 150, 200, 250, 300, 400, 600, and 850 Ha* The IGV total- 
pressure variation was greater than 10% peak-to-peak about the mean total 
pressure for all the frequencies. All the simulation runs were initiated 
from die steady- state operating point on die 100% speed line defined by the 
corrected mass flow rate of 19**^5 Kg/sec (43* 5 Ibj^/sec* ). The settling, 
ramping of P^ amplitude, and steady operation of the model were performed 
according to the model operation procedure outlined in Subsection 4* 1* 1* 

The results of the analysis for the P^G frequencies of 100, 250, and 
600 Hz are summarized »n Figures 22 through 30* Figures 22 through 27 
shew pressure waveforms at various locations of the J85-13/P^G test 
installation. The locations are defined in Table 2* These figures also 
contain the plots of the physical mass flow rate at the IGV and the overall 
pressure ratio across the compressor as a function of time. Figures 28 
through 30 show die blade lift force as a funcrion of time. Similar plots for 
other frequencies, namely 150, 200, 300, 400, and 850 Hz, are presented in 
Appendix C. 

An interesting observation is that the model operation is relatively un- 
stable for the frequencies ranging from 200 to 300 Hz. These frequencies 
are close to the resonance frequency, approximately 230 Hz, of the duct 
between the P^G and the first rotor of the compressor. Therefore, it is 
believed, any small perturbations due to numerical inaccuracy amplify and 
cause a slightly nonstationary operation of the model. When this type of 
nonstationary behavior occurs, the operating points of the compressor cannot 
attain a constant orbit such as the one depicted in Figure 21. Instead, the 
orbit shifts slightly as the average operating point fluctuates sinusoidally 
about a mean value. The frequency of this fluctuation could not be 
determined due to the insufficient amount of computed data. However, this 
type of oscillatory behavior is not severe enough to prohibit operation of 
the model. The simulation can be continued and run indefinitely without 
encountering a major shift of the orbit or drifting into compressor instability* 


61 







I’SACl, KN/m*' Pfiuv, KN/tn" PTAPl, KN/m*’ PTFP6, KW/« 







0 0 012 0 024 


rise, seconds 


cu 



0 0 012 0 024 


rine, seconds 


(U 



Time, seconds 



Time, seconds 




64 


* 

- 1 

. 

m 

m 

S. — 



m 

m 

VT \J 


0 012 0 024 


Tine 9 seconds 


( 

;-| 

JTS 

L 

~i 1 

■1 1 

1 

^ 

! 

i 

• 1 

^ ^ 
1 ■ 

— i 



! 



0 012 0 024 


Tine , seconds 


\ 

- . j 

f 

i 

^ ^ ^ 

1 

^ ! 

1 -J 

1 

^ 

— — 



0 012 0 024 

Tine, seconds 


7\ 

Ay 

VaI A , 

J 

\ R 

JJ 

J V 

M 

m 

V > 

' ! 
t 


0 012 0 024 

Time, seconds 




fp 3 ^ = 250 Hz. 


9 






PSACl^ KN/m"* PTIGV, KN/m PTAPl, KN/m PirP6| KN/« 


cu 


v» — 

^ s_ 


HI 

m 

M 

: V y V 

■ j 

i 

H 

a| 


0 

fU^ 

0* 

Tiott 

— 

012 0 

, seconds 

024 

E 

i 

2! 

H 

OU m 

aAa 

HR 

m 

A/ 

/ viV \ 

1 

1 

1 - ~ 

jgj 

m 

J v 

0 0 
Time 

C\J 

— 

012 0 

, seconds 

024 

CM h- 

E 

i 

§ 

C3 _ 

wAA/' 

BBw 

1 

, a! 

V/ y V 

i 

■■ 

V V 

1 

1 

0 

S 

O'- 

0 

Tioic 

012 0 
seconds 

024 

n 

5 

< 

jp in 

Oi iC 

T 

;■ 1 

/k 

'■■i 

nms 



V 


m 

r \r 

rn 

0 

0 

012 0 

' 024 


Time , seconds 



N m 
< 

CM 

S! 

2 







1 




^ ' 
0 

o -O- 

1 1 

0 

Tim 

[ 1 

012 0 

e , seconds 

' 024 

R 

9 

• 

A A A 

- ! 

♦ 

i 

A A A 

1 A 

7\/V^ 

i 

I ' 

yVV] 

1 

J 

f 

- J1-: 

0 0 

012 0 

' 024 


Time , seconds 


-rr-1 

o/V 




j 





T 1 

0 0 

'012 0 

' 024 


Time I seconds 


= 250Hz (Concluded). 


Figure 25* J85-13 Simulation Results, 


65 



PSrP5i KN/w PSFP3, KN/m PSFP1» KN/m*' P8EMB| KN/m 









P8AC1, KN/m PTIGV, KH/m^ PTAPlt KN/m*" PTPP6, KN/ni 














n 01 m ^ mvo r-co 


UlmUUUUUU 


oqoooopQ 


2 

• I I I I i • I 


OOA^KO«« 






The resonance behavior of the duct betv *' .‘r the P^G and the first rotor 
is shown in Figure ?1. In this figure, file amplitude of the IGV total*pressure 
oscillation is plotted as a function of file P^G frequency The total-pressure 
loss and the amplitude of the P^G area fluctuation are d constant. EIxami- 
nation of this figure reveals that while the same amplitude perturbations are 
introduced at the P^G, the maximum amplitude of the total-pressure 
oscillation occurs at about 230 Ha. Tliis frequency is conunensurate wi*’h the 
resonance frequency computed according to file formula for a closed-closed 
organ pipe. It should be noted that a rotor generates a substantial change in 
acoustic impedence and thus, can be treated as a closed end in an acoustic 
analysis. The P^G plane is choked at all times. Therefore, the duct between 
the P^G and the first rotor can be treated essentially as a closed-closed 
system. 

After the completion of the time-dependent analysis, the model was 
throttled to the point of instability for the purpose of estimating the stability 
}'mit of the compressor under planar inlet distortions. This study is des- 
cribed in the following subsection. 


4. 1. 4 Surge Line Estinnation 


The stability limits of t' e "Mehalic" configuration of the J85 compres- 
sor were obtained by throttling the model to instability (surge). The surge 
line estimation stud^ was conducted at 100% corrected speed for five 
different frequencies, i. e. , 100, 150, 200, 250, and 300 Hz. 


The instability point of the model was determined through observations 
of the migration of the mean operating point. This migration charted on the 
steady-state compressor map was normally well behaved until the point of 
mode] instability was encountered. At the instabil'. point, however, the 
average operating point departed abruptly from the oath it followed and 
wandered erratically until the solution failed. Onc»: the model reached the 
system stability limit, this erratic behavior cf the mean operating point was 
observed even if the throttling was terminated. Therefore, the instability 
point of the model was determined as the point ct which the average operating 
point abruptly departed from its path even in the absence of throttling. 

After the point of instability was defined, the system model was 
throttled to instability to investigate the effect of the P^G frequency. First, 
the model was brought to steady operation through settling and ramping after 
it entered computation at s 19.79 Kg /sec (43.5 Ib^/s^c). The amplitude 
of the P^G area oscillation was 10% peak-to-peak for all frequencies. The 
total -pressure loss across the P^G was set at 45% of the inlet total pressure. 
The loss is commensurate with fiie actual total -pres sure loss across the 
P^G determined experimentally by NASA. After the mod^l was allowed to 
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run for approximately 10 to 20 dycles of the planar wave, the model was 
throttled for a period of 2 cycles. Hien, die model was run for 4 cycles and 
observed for occurrence of instability. If die model was stable, it was 
dirottled again for 2 cycles and then, observed for 4 cycles. This process 
was repeated until die point of model instability was observed. The throttling 
was conducted at a rate slow enough to assure a quasi steady«state dirottUng 
process; die exit flow function was reduced at the rate of 0. 00005 (^^-ft) 
or less. 

The results oi the stability analyses are presented graphically in 
Figures 32 throu<>a 36. These figures show die instability point, as well as 
the locii of die mean operating point for ramping and throttling. EIxamination 
of these figures reveal that the locus of die average operating point deviates 
furdier from the steady<>state speed line at P-^G Trequencies that are close to 
the resonance frequency of die P^-to*first compressor rotor system. It is 
believed that the larger deviation is due to the large amplitude of the planar 
wave induced by die resonance of the duct (see Figure 31). 

The loss in the ^mrge pressure ratio is also greater at diese frequen- 
cies. This is illustrated clearly in Figures 37 and 38. These figures present 
the distortion sensitivity as a function of the P^G frequency and the nor- 
malized frequency, f*, defined as 

_ frequency of P^G 

resonance frequency of P^G-to-first rotor system 

The distortion sensitivity is defined by the following relationship 

Ds = Aprs/ Ap^/zp^) 

where APRS and represent the loss of surge pressure ratio and the 
average total pressure, respectively. Figures 37 and 38 show that the distor> 
tioa sensitivity is maximum at approximately fp3Q = 230 Hz or f* = 1. 

4. 2 FREX3UENCY RESPONSE ANALYSIS 


The frequency response characteristics of the J85-13/P^G tesi 
installation were investigated using the Digital Surge Prediction (DSP) pro- 
gram described previously. The model geometry employed in this linearized 
analysis was essentially the same as the one used in the time -dependent 
analysis. However, the geometry was broken into two parts, i. e. , 
bellmouth-to-P^G and P^G-to-combustor exit. The frequency response 
analys's was performed separately for the two sections. Further, the 
number of volumes were -educed to 28 and 31 for the sections ahead and aft 
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of the P^G, respectively (59 total); v^ereas 79 volumes were used to describe 
the entire system in tiie time«dependent analysis. These modifications to the 
geometry were necessary because of the limited memory capacity of the 
available computer and exponential overflow and underflow in evaluating file 
characteristic equation of a large size matrix. 

For file section ahead of file P^, fiie analysis was conducted for the 
frequencies ranging from 10 to 600 Hz. The average lengfii of the volumes 
describing this section was 4 inches. Therefore, 600 Hz was the maximum 
frequency at ^ich the model can be expected to yield reasonably accurate 
response characteristics. The frequency range for the region aft of file P^G 
was from 10 to 1000 Hz. 

The boundary conditions for the bellmouth>to-P^ system were: 1) 
Constant total pressure and total temperature at file bellmouth, and Z) 
Oscillating flow function at the P^G. In an actual engine test, the entire 
test setup or the bellmouth is placed in a large plenum. Usually, the 
total pressure and total temperature of the plenum are not affected by the 
operation of the ei^ine. Therefore, the constant total-property condition was 
employed as the inlet boundary conditions. The P^G is considered as a device 
that primarily oscillates the mass flow rate. However, the DSP program 
does not allow fluctuating the exit flow. Thus, an approximation of the 
actual boundary condition, i. e. oscillatit^; flow function, was introduced at 
the exit. The frequency analysis based upon this approximate exit boundary 
condition may yield amplitudes and phase angles which differ from the actual 
values, whereas, it should give correct resonance frequency and locations of 
nodal points. Consequently, the results of the frequency analysis for the 
section ahead of the P^G should be used only in studying the frequency 
response trend. 

The boundary conditions for the P^G-to-combustor exit system were; 

1) Oscillating flow at the P^G, and 2) Constant flow function at the combustor 
exit. As stated previously, the P^G oscillates the mass flow rate. There- 
fore, oscillating flow was introduced as the inlet boundary condition. During 
the actual operation of a gas turbine engine, the turbine diaphragm is choked 
due to high pressure and high temperature of the combustor gas. This 
choked condition was siniulated by assigning the constant flow function to the 
combustor exit plane. These two boimdary conditions correctly describe the 
actual physical conditions at the respective planes. Therefore, the resultant 
frequency analysis for the section aft of the P^G will be exact wi.hin the 
accuracy of computation. 

Prior to performing the analysis at various frequencies, the frequency 
response characteristics for the section aft of the P^G at 100 Hz were 
compared with those obtained from the time -dependent analysis co verify the 
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validity of the DSP program. The time -dependent analysis was performed 
with less tLr.n 0. 5% IGV total-pressure variation about the mean value to 
approximate the small perturbation condition of tiie linearized analysis. The 
comparison shows that ^e static pressure response in tiie duct between &e 
P^G and the IGV of the compressor computed by tiie two different methods 
agrees well in amplitude. Figtire 39* as well as in phase angle. Figure 40. 

The amplitudes calculated by the DSP program differ from diose obtained by 
the time -dependent analysis by less than 10% at all locations. The phase 
angles show a closer agreement. Therefore, it is concluded that the 
linearized analysis technique can be employed successfully in evaluatii^^ the 
resonance frequency, nodal points, etc. of a compression system such as 
the J85-13 /p3g. 

A similar comparison and the resultant excellent agreement between 
the linearized analysis and the time -dependent analysis in evaluating the 
frequency response characteristics of the J85-13 engine are presented in 
Reference 4. Therefore, they are not discussed here. 

Subsequently, the frequency response analysis for the two sections 
were performed at the frequencies stated previously. Some of the results 
are summarized in Figures 41 through 48. The remaining results are 
presented in Appendix D. The results show that the dominant resonance 
frequency of the P^G-to-combustor exit system is approximately 230 Hz, 
Figure 41. This value is equivalent to the resonance frequency of the duct 
between the P^G and the first rotor of the compressor evaluated according 
to the closed-closed formula for organ pipes. Figure 42 shows that the 
bellmouth-to-P^G system has pressure response peaks at about 31, 90, 150, 
210, 250, and 280 Hz. Note that 31 Hz is the fundamental resonance frequency 
of the bellmouth-to-P^G system and 90, 150, 210, and 280 Hz are approxi- 
mately equal to the odd harmonics of the open-closed system. It is believed 
that 250 Hz is associated with another open-closed system within the 
bellmouth-to-P^G section. This subsystem is defined by the geometry 
between the second strut and the P^G, see Figures 1 and 10. 

Figures 43 through 48 present the results of the nodal analysis for the 
two sections of the J85-13/P^G test installation at different frequencies. 

4.3 VERIFICATION STUDY OF THE EXPLICIT ANALYTICAL SURGE 

CRITERIA 


The analytical surge criteria derived in Subsection 3. 5 were tested 
with the "Moss" and "Mehalic" configurations of the J85 engine. The criteria 
yielded instability points slightly beyond the experimental surge line for both 
engines at 100% speed. However, these criteria were ineffective at low 
speeds for both configurations of the engine. It is believed that the 
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ineffectiveness of tlie criteria at low speeds is mostly due to die inadequate 
boundary condition employed at the rotor exit. It was assumed diat die flow 
fiinction at each rotor exit is constant, lliis assumption is adequate at 100% 
spc i since die corrected flow rate does not change significantly along the 
speed line; dius, the relative Mach number of flow at the exit of eadi rotor 
is almost constant. However, diis boundary condition is inappropriate at 
low speeds since die corrected flow rate and consequently, die relative Mach 
number changes significantly as the operating point migrates along a speed 
line. 


Figure 49 presents die instability point of the "Moss" ei^ine at 100% 
speed as computed by die explicit analytical stability criteria along with the 
stability limits determined by other methods. Figure 50 is a similar plot for 
die "Mehalic" engine. Hie occurrence of the instability was indicated by the 
second inequality. Equation 76, of the criteria; the first criterion, Ek)uation 
75, was satisfied even when the second criterion indicated an instability. The 
second criterion contains die slopes of the blade characteristics, a>', 5 and 
fict see Appendix A. It is believed diat these slopes play an important role 
in the stability of the compressor and thus, the instability was indicated more 
readily by the second inequality. For both engines, die instability occurred 
at the first rotor only. 

A critique of the explicit analytical stability criteria is presented 
through a discussion of the three major assumptions employed in their 
derivation. These assumptions are: 

1) The flow is homentropic. 

i) At each rotor exit, the flow function is constant, and 

3) The stability of die entire compressor can be determined by 
examining the stability of individual rotors. 

The first assumption was introduced to simplify the derivation and thus, 
the resultant criteria. This condition allows one to write Equation 69, 
Subsection 3. 5, as a 2 x 2 matrix instead of a 3 x 3 matrix since the energy 
equation vanishes. Also, the characteristic equation. Equation 74, is 
reduced to a quadratic equation from a cubic equation. Consequently, 
derivation of the criteria becomes simpler. It was assumed that the vari- 
ation of the entropy has an insignificant influence on the compressor in- 
stability. However, it is clear that the exclusion of the entropy as a para- 
meter does not allow the criteria to detect any instabilities that may be 
indicated by the violation of the second law of the thermodynamics. 

Introduction of the second assumption as a boundary condition was dis- 
cussed before. It was stated that this assumption is quite valid at high 
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speed; however> it is not nearly so valid at low speed. The inadequacy of 
this boundary condition at low speeds is demonstrated clearly by the fact diat 
no instability point was predicted by the criteria at these speeds. 

The third assumption was introduced also to simplify die explicit cri« 
teria. If the criteria are applied to the entire compressor, the character- 
istic equation of a large size matrix should be evaluated and consequently, 
derivation of explicit stability criteria will be extremely difficult, if not 
impossible. The stability criteria for each rotor volume, however, cannot 
account for the interaction of the instability of other rotor volumes widi that 
of the rotor volume under consideration. 

Hie feasibility of employing explicit analytical stability criteria to 
predict the instability was demonstrated. This was accomplished by obtain- 
ing the instability points of two different configurations of die J85 engine at 
100% speed that agree well with the experimentally determined stability 
limits. However, it is believed that a furdier investigation, and ensiing 
improvements, of the analytical stability criteria is necessary to obtain more 
reliable stability predictions at all operating conditions of a compressor* 

This investigation should include: 1) Introduction of the energy equation in 
the form of entropy change into the derivation of the criteria to account for 
the possible violation of the second law of the thermodynamics at the insta- 
bility limits, and 2) Development of boundary conditions or other techniques 
that allow for a stronger interaction of compressor stages in evaluating 
stability. In addition, it is desirable to explore the possibility of obtaining 
stability criteria through a completely new approach such as a parametric 
study that relates the occurrence of instability to the geometric and aero- 
dynamic parameters of a compressor. 
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5.0 CONCLUSIONS AND RECOMMENDATIONS 


The main objective of this program was to develop a computer simu* 
lation program that is capable of modelling the operation of the J85-13 
Planar Pressure Pulse Generator test installation. This task was accom- 
plished successfully by modifying an existii^ computer program tiiat was 
developed to simulate a test installaticn» called die FI 01 Fan/P^G» similar 
to die J85-13/P^G. This conclusion is supported by the following specific 
a .complishments: 

1) An expression for the unsteady-lift force was developed and added 
to the existing program to improve the modelling of die unsteady 
aerodynamic behavior of compressor blades, 

2) The modified program was validated dirough a comparison of the 
simulation results with the experimental data obtained during the 
FlOl Fan/P^G test, 

3) The simulation program was utilized successfully in modelling 
operating point excursions of the J85-13 engine about a stationary 
mean operating point. The model was s^*lble over a wide range of 
P^G frequencies and amplitudes, and 

4) The J85-13 engine was then throttled to the stability limits and the 
planar wave distortion sensitivity of the engine was determined at 
various P^G frequencies. The instability points were identified by 
a consistent pattern of the model behavior, i. e. , a rapid departure 
of the mean operating point from a well-behaved path. 

As supported by the items listed in the preceding paragraph, the 
J85-13/P^G simulation program was developed successfully and its validity 
was verified. However, due to lack of proper interstage and over -the -rotor 
test data, the validation of using the unsteady-lift force for the description 
of the unsteady aerodynamic behavior of compressor blades was indirect, 
i.e. , through a comparison of the analytical and experimental pressure 
waveforms in the ducts ahead and aft of the FI 01 fan. Therefore, it is 
recommended that these types of test data be taken during the future 
J85-13/P^G test and be utilized for the direct validation of the simulation 
technique. Such experimental data will be useful also in refining the expres- 
sion for the unsteady-blade-lift force as well as in evaluating the constants 
contained in the expression. 

The frequency response characteristics of the J85-13/P^G test 
apparatus were determined in support of selecting instrumentation locations. 
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The fiindamental :resonance frequencies of two subsystems of the test instal* 
lation, i. e. , bellmou^-to-P^G and P^G«to>combustor exit, were identified. 
The nodal analyses of tite two subsystems were performed as well. 


In addition, the feasibility of employing a set of explicit analytical 
criteria in evaluating the stability of an axial compressor was investigated. 
These criteria successfully yielded the instability points for two different 
configurations, "Moss" and "Mehalic", of the J85-13 compressor at 100% 
speed. However, they were ineffective at lower speeds. It is believed that 
this was mainly due to the boundary conditions whose validity was weakened 
at lower speeds. The following items are recommended to improve the 
explicit analytical stability criteria. 

1) Development of new boundary conditions that describe the acual 
aerodynamics more accurately at low speeds, 

2) Development of a new technique that permits a proper coupling of 
stages in determining the instability, 

3) Addition of entropy as a parameter to account for the possible 
violation of the second law of thermodynamics of instability 
points, and 

4) Feasibility study of alternative approaches for the derivation of the 
stability criteria. One of such approaches is a parametric study 
that relates the occurrence of instability to the geometric and 
aerodynamic characteristics of a compressor. This approach may 
yield an independent set of criteria as well as serve as a means of 
enhancing and simplifying the existing criteria. 
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* APPENDIX A 


Deriv^on cf Explicit Analytical Surge Criteria 
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In this appendix, the relationships used in the explicit analytical 
criteria (see Sub-section 3. 4) are derived. 

The governing equations for the homentropic flow of a compressible 
fluid through a compression system are given in volume-averaged form as. 
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For this nonlinear system, the stability criteria are given as 
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where the partial derivatives are 
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The partial derivatives in the right>hand side of E2quations A-7 through 10 
are derived below. 
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The relation between pressure and density for an isentropic process 
of a perfect gas is 
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where C is a constant. Thus, the partial derivative, can be written 

as 1 
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The nnass flow rate, can be expressed in terms of density and 

velocity as 
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where A, is the area and C « is the axial velocity. 
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For an i sent topic flow, the density ratio is 


1 



Thus, the derivative din becomes 
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Substitution of Elquation A>19 and the i sent topic relation 
dlnP, 
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into Equation A- 17 yields 
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Subsequent introduction of this relationship into Ejquation A- 15 gives 
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The boundary condition at the exit of a rotor volume is given as 
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where C is a constant. 
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Substitution of the isentropic relations 


P, =CAT 




into Equation A>23 and rearrangemeoft yields 
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Pj > cw^ 


y* 1 


Therefore, the partial derivative becomes 
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Equation A>26 can be rewritten, with the aid of the relation given by 
Equation A -24. as 


P^ = C W^ 
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Therefore, the partial derivation of P^ with respect to is 
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Derivation of 
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The blade force, Fz, can be expressed in terms of the tangential 
force, F-y, and axial component of the drag force, 


F* = Ft 



(A- 30) 
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where the corrected>Uft •direction angle* and the drag force* F_ * 

are not a fimctioa of die exit independent variable* W 2 . Therefore* ^ 
the partial differentiation of widi respect to W 2 yields 


dF^ dFj 

avf^ ~ aw 2 


(A-31) 


From Euler's Turbine Equation* the tangential force is obtained as 

rZ 


- I r ^2*1 ijn, 1 
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where r j and are the pitch-line radii at the inlet and exit of a rotor, 
resoectivelv. 



Figure A-2. Rotor Exit Velocity Diagram^ 

Introduction of the angular relationship, depicted in Figure A-2, 

“z^z 

tano, = ^ - tan/J, (A-33) 

c w 2 2 

enables one to rewrite Equation A-32 as 
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The third term in the right ‘hand side of Equation A-34 is independent of 
Further, the relative eadt air angle, and the pitch-line blade speed. W^, 
are also independent of Therefore, the partial derivative, 

becomes 
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Substitution of this relationship transforms Equation A-31 into 
dF 2r , W 
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All the quar4tities in the right-hand side of Equation A-30 are a function 
of the density, Pj. Therefore, 
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Derivation of the partial derivatives in Equation A-37 is lengthy and 
tedious. Therefore, only the results are presented here. 
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The slopes, tn^ and of blade characteristics are evaluated 

from the polynomial representatives of the loss coefficient, cu', the deviation 
angle, tan 5, and the lift -direction correction angle, tan/3^. 'Diese poly- 
nomials as a function of tani are discussed in Section 3. 2. 
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APPENDIX B 

Results of FI 01 Fan/P^G Sirrmlation 


109 







Figure B-2. FlOl-Fan Blade Row Total-Pressure Ratio as a 
of Time, 100% NA/e*, fp3^ = 80Hz. 








Time, seconds 


Fioure B^4. FIOl-Fan Blade Row Total -Pres sure Ratio as a Function 
of Time, 100% N/JT, fp3^ = 220Hz. 
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Figure &-6. FlOl<*Fan Blade Row Total-Deflq>erature Ratio as a Fiinction 
of Tine, 100% N,'/r, fp3^ . 42Hz. 
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rise , seconds 


Figure B-8* FlOl-Fan Blade Row Total*Teaq>erature Ratio as a Function 
of Ti«e, 100» N/V*", fp3g * 118 Kz. 
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Figure E-9« FlOl^Fan Blade Row Total -Temperature Ratio as a Function 
of Time, 100% H/JT^ fp3g = 220Hz. 
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. PlOl^Fan Blade Row Total-Teaperature Ratio as a Function 
of Ti«e, 100% fp3g = 350 lfa. 
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APPENDIX C 

Results of J85-13/P^G Simulation 
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Figure C-1. J85*13/P^G Simulati 
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Figure C-3. J85-13/P^G 
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Figure C-4. J85-13/P^G Simulation Results, = 200Hz (Concluded).. 
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Figure C-8. J85-13/P^G Simulation Results^ (Concluded) 
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Figure C-9. J85-13/P^G Simulation 
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APPENDIX D 

Results of J85-13/P^ Frequency Response Analysis 
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